Reheating After Quintessential Inflation and Gravitational Waves 
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We investigate the dependence of the gravitational wave spectrum from quintessen- 
tial inflation on the reheating process. We consider two extreme reheating processes. 
One is the gravitational reheating by particle creation in the expanding universe in 
which the beginning of the radiation dominated epoch is delayed due to the presence 
of the epoch of domination of the kinetic energy of the inflaton (kination). The other 
is the instant preheating considered by Felder et al. in which the Universe becomes 
radiation dominated soon after the end of inflation. We find that the spectrum 
of the gravitational waves at ~ 100 MHz is quite sensitive to the reheating process. 
This result is not limited to quintessential inflation but applicable to various inflation 
models. Conversely, the detection or non-detection of primordial gravitational waves 
at '^lOO MHz would provide useful information regarding the reheating process in 
inflation. 
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I. INTRODUCTION 

Inflationary paradigm solves various problems in the standard big-bang model (flatness, 
horizon, monopole problems, etc.) and predicts density perturbation and gravitational wave 
spectra which are both almost scale-invariant. The two of three major predictions of in- 
flation (flat universe, scale-invariant density perturbations) are consistent with the current 
observations, and the detection of primordial gravitational waves is a litmus test of the in- 
flationary paradigm. The dynamics of inflation is believed to be described by a scalar field 
called "inflaton" (or multiple of it) with the energy scale of ~ lO^^GeV. 

On the other hand, recent observations indicate that the present Universe is dominated 
by dark energy: an (approximately) smooth energy component with negative pressure. If 
dark energy is dynamical, its dynamics is likely to be described by an ultra-light scalar field 
with mass ^ Hq, called "quintessence". 

Then, it is natural to imagine whether it is possible to unify the inflaton and quintessential 
fields. A toy model was proposed by Peebles and Vilenkin in which the potential consists of 
two parts: + (0 < 0) for inflation and AM7(0^ + M^) (0 > 0) for quintessence Q. 
The two potentials are glued at = by hand. In this paper, although not essential, as a 
slightly elaborated model, we consider a smooth, exponential type potential that describes 
both the inflaton and quintessential parts^ in the hope that the model could be embedded 
in particle physics framework. For other models of quintessential inflation, see P]. 

In order to recover the hot big-bang universe, the universe has to be reheated after infla- 
tion. In models of quintessential inflation proposed so far, the only mechanism of reheating 
that has been considered is the gravitational particle production in an expanding universe. 
The spectrum of gravitational waves based on the gravitational reheating was calculated 
js, y, 0] and a strong enhancement at ~GIIz was found (see also in different context). 
However, the gravitational reheating process is very inefficient and may lead to cosmological 
problems associated with overproduction of dangerous relics such as gravitinos and mod- 
uli 0]. Recently, Felder, Kofman and Linde proposed a much more efficient mechanism 

paper, we calculate the spectra of 
gravitational waves in both cases of the instant preheating and gravitational reheating, and 



^ This possibility was noted in a footnote of P|. 
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examine the dependence of the gravitational wave spectrum on the reheating process. 

This paper is organized as follows. In Sec. II, we introduce our model of quintessential 
inflation. In Sec. Ill, we briefly review two reheating processes after quintessential inflation. 
In Sec. IV, we calculate the spectrum of the gravitational waves. Sec. V is devoted to 
conclusion. We use the units of ^ = c = 1. 



II. YET ANOTHER QUINTESSENTIAL INFLATION 

We utilize an exponential potential of the form (Mpi = 1/ VSnG = 2 x 10^^ GeV) 

\/(0) = K)exp(-A0/Mpi), (1) 

but with A being a function of 0; A = A(0). 

Inflation. If the Universe is dominated by a scalar field, then the scale factor a(t) be- 
comes a oc t^/^^ and the power-law inflation is realized for A < y/2 j^. Recent observations 
of cosmic microwave background anisotropies indicate that the spectral index of the scalar 
perturbation, n<j, is consistent with flat one: = 0.93±0.03 'l^. Since n^ — 1 = — A^ under 
the slow-roll approximation, we require 

A < 0.3. (2) 

If we assume that inflation terminates at (p Ohj changing the slope A, then the amplitude 
of the density perturbations observed by COBE fixes the energy scale of inflation: Vq ~ 
10-i2MpV 

Kination. After the end of inflation, the universe is still dominated by 0. Inflation 
is followed by a kinetic energy-dominated epoch (hereafter called kination). During the 
kination, 0^ ~ Vo{a/akd)~^ and the scalar field evolves as 

~ V6Mpiln(a/afcrf), (3) 

where akd is the scale factor at the end of inflation. Therefore in order for the potential energy 
V oc exp(— A0/Mpi) not to dominate the energy density of the Universe, it is required that 

X>VQ. (4) 

The presence of the epoch of kinetic energy domination is the main reason of enhancement 
of amplitude of gravitational waves in quintessential inflation models discussed below Q . 
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Dark Energy Domination. The Universe must eventually be dominated again by the 
scalar field to account for the acceleration of the Universe. This is accomplished by the 
termination of the domination of the kinetic energy of the scalar field and hence we require 

A < ^6. (5) 

The scalar field stops moving when V{(f)x) ~ -^o^pi the Universe becomes eventually 
dominated by the scalar field again. 

From the requirements (H)) and (0), we find that different three slopes (A) are needed 
for the stages of inflation, kination and dark-energy dominance. As an example, we consider 
a potential in the following form: 

\/(0) = K)exp(-A(0)0/Mpi), (6) 

where A,u and v are constants. In the following analysis, we take A = 1.1 and arrange u 
and V to flx Q^p = 0.7 today. Inflation is terminated when {aH)~^ is minimum. Vq is flxed 
by the amplitude of density perturbation observed at COBE scale : V^^'^ / {'\/75tiM^]V') = 
1.9 X 10"^ 

Possible Scenario. We speculate on possible scenarios for realizing the above potential 
form. The flrst one is the system of the Lagrangian consisting of a non-canonical and a 
simple exponential potential: 

- iz($)2(V<l>)' - Voexp{-a^). (8) 



A non-canonical kinetic term appears in supergravity theories [Ul and in higher-dimensional 

theories via dimensional reduction. Introducing a new variable such that = 

so that = / = we may obtain the required form of the potential with 

n n 

Another scenario is a scalar fleld with non-minimal coupling with gravity |12l . Il3| : 

f{<P)R-l{V<j>r-V{<f>). (9) 
Since the effective potential for the dynamics of the scalar fleld is V{(f))/ /"^{(p) rather than 



12l llS^ . we may realize the required form from a simple form of V{(f)), although we 
must be careful about possible time-variation of the gravitational constant and deviations 
from general relativity 1^ . 
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III. REHEATING PROCESSES 

Gravitational particle production. Massless particles are created due to the change of the 
expansion-law of the space-time. These particles can reheat the universe. Because the energy 
of massless particles decreases like that of radiation, the energy density of created particles 
eventually dominates even if it is initially negligible. The created particles thermalize the 
Universe in the following standard process [l]. The energy density of the created particles 
is pm = RH^di^/^kd)~^, where H^d is the Hubble parameter at the end of inflation and 
R ~ 10"^ per scalar field component 3] and R ~ 10~^A^s for A^^ scalar fields. 

The thermalization of pm occurs when the interaction rate na is faster than the expansion 
rate of the Universe H. Here the number density n and the cross section a are written in 
terms of the energy of created particle e ~ Hkdid/cLkd)'^ as ~ Re^ and a ~ a^e~^ with 
a ~ 10~^. Thus the thermalization occurs at ath/dkd ~ a^^R^^^'^ and the temperature is 
Ttr, = R^/^aHkd. 

Since pmj P4> ~ RHldM'i^ {a / akd)'"^ during kination, the radiation dominated epoch be- 
gins at 

<^rlakd = ^|lHkdM,, ~ W'R-y^ {j^i^) ■ (10) 
From this result, we can obtain the temperature at the beginning of the radiation. 

Tr ~ i?3/4^^^(3^2 )-l/2 _ ^q6^3/4 f J^M \ ^^-^ ^^^^ 



VlOi2 GeV. 

This temperature is sufficient for the beginning of the standard hot big bang model. A 
typical result of the time evolution is given in Fig. 1. 

Instant preheating. We consider the so-called instant preheating as another reheating 
mechanism. Felder et al. proposed a much more efficient mechanism of reheating which 

n 

works even for a potential^ without a minimum 7]. The mechanism is similar to preheating 
by parametric resonance [8]. They introduced the interaction Lagrangian, —\g'^ (f^X^ '~\'4''4'Xi 
where g and h are the coupling constants. After the end of inflation, the energy of the inflaton 
is transferred to the scalar field x and the field x decays to the fermions ip. 

For simplicity, we assume that infiation ends at = and the scalar field x does not 
have a bare mass. Then, the effective mass of the field x is f^ix ~ 9\'P\- The particle creation 
of the field x occurs when the adiabaticity condition is violated, > m^. This happens 
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right after the kination begins. Hence when 



\<l>M\/9^V,'^y^, (12) 



where <pkd is at the end of inflation. In order for the particle creation to be completed 
within a sufficiently short period, we must have < Mpi and g must satisfy g > 10^^. The 
time interval of the particle production is then written as 



At ~ -p- - g-'^X-'^'. (13) 
\(pkd\ 

Thus the created particles will have typical momenta k ~ At~^ ~ g^^'^Vo^^ 0- The oc- 
cupation number of the created particles is n{k) = exp (^—'^'k^/g^/\^^ j^. Integrating this 
equation, we can get the number density of the field x 

f — \ 

n^ = ^ dkk^n{k) = ^ / . (14) 



27r2 Jo ' ' 87r3 

The energy density immediately after the creation is p^^ ~ (s'xA^) /Svr^. Soon after that, 
the field x decays into the fermion ip with decay rate F = h'^m^/Sir = h'^g\(f)\/8'JT. The 
subsequent decay of the field x leads to a complete reheating of the Universe. 
The equation of the field (p after particle production is written as 

+ 3i/0 + /(x')0 = O, (15) 

where the potential can be neglected in the kination. When the field becomes greater 
than 0^ and particle production ends, the particles x become nonrelativistic. Then we can 
calculate (x^) 

]_ f rikk'^dk _ % _ {gVVo)^^^ fakdV 



' 27r2 ] .JWTgW g<P 87r3^0 \ a J ^ ^ 

Therefore the backreaction of created particles on is negligible if g^ (x^) < 3H(p, that is 

«<«^^|^^I7^j • (17) 

We consider the case of ab/akd > 1, so we get g < 0.8 from Eq. (fT7j) . If F > Hi,, where Hf, is 
the Hubble parameter at the scale factor ai, = a^d [8VQTTX^y{g^/^Mpi)Y/^, particles x will 
decay to fermions ip a.t a < at, and the force driving the field back to = will disappear. 
This condition is given by 

> (18) 
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FIG. 1: Evolution of the energy densities for gravitational reheating. The solid line is the potential 
energy density of (f). The gray line is the kinetic energy density of (p. The dotted line is the radiation 
and dust energy densities. 
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FIG. 2: Evolution of the energy densities for instant preheating. The solid line is the kinetic energy 
density of (f), the dotted line is the energy density and the dashed line is energy density of ip. In 
the right panel there is the era dominated by the scalar field X: because the decay rate of the field 
X is low. 



where we have assumed (f) ~ Mpi. When the backreaction is nonnegligible, the scalar field 
X dominated epoch begins. The reason is that the condition (|17p is essentially the same as 
that for < p^. When if ~ F, the field x decays into ip and the energy is transferred to 
the radiation energy. Thus the radiation dominated epoch starts. In Fig. 2, we show the 
results of time evolution of fields. 
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IV. GRAVITATIONAL WAVES 

In each reheating process, we take A = 1.1 in Eq.© and arrange u, v to fix = 0.7 at 
the present. From this result, we obtain the scale factor and the Hubble parameter at the 
beginning of each epoch. 

We calculate the graviton spectrum produced by an infiationary stage by using the Bo- 



goliubov transformation in terms of conformal time r [17|, Q| . But the spectrum shape is 
easily understood in the following way before carrying out a detailed calculation. When 
the gravitational waves come out of the horizon during infiation, the amplitude of the grav- 
itational waves is frozen and the value is proportional to Hout, where Hout is the Hubble 
parameter when the scale of the perturbation crosses the Hubble horizon. After reentering 
the horizon at scale factor aj„, the amplitude decreases as 1/a. So the today characteristic 
amplitude hc{uj) which reentered the horizon in an era with scale factor a oc is 



^ ^ fi3/2)-u+in/n-l) 

'"c ^ ^^oui 

where 



/.,ociJ„„,^oc/(=^/2)-'^+("/"-^). (19) 
ao 



^'l + TT^)^ (20) 



2 V* ' 2- 

and A here is its value during infiation which we approximate by A(0 — oo) in Eq. (jH)) 
for simplicity. This means that the scale factor during infiation is approximated to be 
proportional to t^/^^ . From this, the energy density of the gravitational waves per unit 
logarithmic interval of frequency is 

where pc is the critical energy density. Our model consists of the five phases: infiationary 
phase (IN), kinetic-energy dominated phase (KD), radiation dominated phase (RD), matter 
dominated phase (MD) and potential dominated phase (PD). Among them, the phases 
when the gravitational waves reenter the horizon are KD (n = 1/3), RD {n = 1/2) and MD 
{n = 2/3). Correspondingly, the spectrum has three branches which are proportional to 
j!i+{3~2iy)^ f3-2u ^j^^ j~2+{3-2u) ^^^^ ^^ie shapes are showed in Fig. 3. 

The detailed calculation using the Bogoliubov transformation is as follows. The scale 
factor in each phase is given by 

/ ^ X (l-2!/)/2 

a{T) =akd[ — for IN, 

\Tkdy 
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a{r) = akd 



:i-2z/)( — + -— 



aid 1 - 2z/ x, / r 



Tkd 



a(r) = Qd 



2 Tkd \rr flfcrf 1-2U Tr 



1 - aid Td 



arfflr T fed 



Td 1 - 2z/ Td 



for KD, 



for RD, 



for MD, 



a r 



1-^ 



1-^ 



1 



forPD, (22) 



2 ^/a^ — yorf 

where subscripts fee?, r, d and p mean the epochs of the phase transitions between IN-KD, 
KD-RD, RD-MD and MD-PD stages, respectively. In the PD phase, we have assumed that 
the universe accelerates exponentially, since the potential energy density is nearly constant 
as seen from Fig. ^ The conformal time at each phase transition is expressed as 

\-2v 1 



O'kdH, 



kd 



a 



kd 



Td 



Od 



2akdHkd 
1 



+ Tkd, 



Gjrp Hip 



ad , 



adH, 



+ rd. 



(23) 



Note that Tkd < 0. 

In the TT (transverse-traceless) gauge, a positive frequency function for the gravitational 
perturbation with comoving wave number k is represented by 



/.;f) = a(r)M?(k)/i,(r)e^'^-, 
where e\j\k.) {cr = 1,2) are two independent polarization tensors normalized as 

e('^)'^'(k)eSf' (k) = 5'^'^'. 



(24) 



(25) 



The Einstein equations give an equation for fik 



IJk" + 2— /ifc' + k'^jjk = 0, 



(26) 



where the prime denotes the derivative with respect to the conformal time. The mode 
function fi satisfies the normalization 



(27) 
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In the inflationary phase, assuming the standard Bunch-Davies vacuum, the solution is given 

by 

/^.(r) = /^(-r,,)^-''(-r)'^iJ«(-fcr), (28) 

where the subscript i means it is the solution in the inflation phase and if^^-* {x) is the Hankel 
function of the flrst kind. 

The solution in the kination is 



VkdiV = ( —Tkd 



H^^^ ( k 



2v 



T 



2v 



Tkd 



(29) 
(30) 



2aM \\-2v 

where a^d and /J^d are called the Bogoliubov coefficients, il}^^ is the mode function in the 
kination normalized as Eq. (|27p. and Hq^\x) is the Hankel function of the second kind. The 
Bogoliubov coefficients are determined by the condition that /i and /i' be continuous at the 
phase transition. From this condition at r = Tkd we get 



Pi 



Hi'\-kna) 

ivr / 4 



1 



1 - 2z/ 



1/2 



(1) 



1 - 2u 



-kri 



kd 



kd 



1 - 2z/ 
k 



1/2 



2 ^ 



(2) 



Tkd 
1 



1 - 2z/ 



kTkd) H^^H-knd) - iyn}HP(-kn,)HP 



'kd 



kd 



ikdjJ^Q 



2v 



-kri 



kd 



+ 2 ° 



(2) 



1 



kud ) iHl'\{-knd) - Hltlii-kTi 



(1) 



kd, 



(31) 



,1 - 2u 

For wavelengths of cosmological interest, we have —krkd 1- When we assume —kTkd ^ 1, 
we get 



Viv) (\-2v\^l^ 



4 



where we have replaced r^d by its expression in terms of the scale factor and the Hubble 
parameter given in Eg. (1231) . 

nn 

Here we need to mention the adiabatic theorem |17l . I18I |. The transition of the phases 
is not instantaneous in reality. The change of the cosmological evolution occurred with the 
time scale AT = . The time scale associated with the mode function with the comoving 
wavenumber k is the frequency uj = k/a. If k/a > AT = H~^, this change of phase is 
adiabatic and no particle creation occurs. This results in a = 1, /3 = for A; > aH at an 
epoch of a phase transition. 
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In order to obtain the spectrum of the present gravitational waves we need the solution 
in the potential dominated epoch. It is acquired by the the Bogoliubov transformation. The 
relation between the mode functions are written using the Bogoliubov coefficients. 

ipd ^ apipp + ppip*. (33) 

These Bogoliubov coefficients are obtained by requiring the continuity condition at each 
phase transition. The mode function at the PD epoch is expressed as 



l^p = Oltotal i>p + Ptotal i>p, 



(34) 



where atotai and jStotai are given by 
/„, ^ \ / 



Oltotal Ptotal 
\ Ptotal 



''total 



Oikd Pi 



\ Pkd ^kd 




Oip Pp 

yl3p Oi*pj 



(35) 



In the PD epoch, the wavelength cannot come into the horizon because of the accelerated 
expansion. So, for all scales which we can observe {k > GoHq), we have ap — 1 and /3p = 0. 
The number density of created gravitons with conformal wave number k is 



Nik) = \Ptotal\ 



(36) 



The energy density of gravitational waves with frequency / = k/2Trao is written taking into 
account of two polarization of the gravitational waves as 



dpgw(/) = 2 (27r/) |Aotair47rfrf/. 



(37) 



Prom this the density parameter of gravitational waves is 



Ogw(/) = 167rVe-'|Aotair/' 



(38) 



Thus we finally obtain 

4(1 - 2iy)r^{u) 



^gw(/) 



2Hkd o-kd 
7r(l — 2y) oo 



2v 



-2v) 



Hr O-r <- y -^fc(i (^kd ^gg^ 

27r Oo 27r oq 
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FIG. 3: The spectrum of gravitational waves. The left panel is in the case of the reheating 
by gravitational particle production. The right one is in the case of instant preheating. For 
comparison, we also plot the result of gravitational reheating with R = 1. 
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'2'Hkd cikd 
7r(l — 2z/) ao 

"^Hkd Clkd 

7r(l — 2z/) 



2u 



Hd CLd „ Hr CLr 

ZTT ao /vr ao 



2u 



ao J ao 
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(41) 



Note that u = ?)/2 corresponds to de Sitter inflation. 

The spectrum of gravitational waves for each reheating process is shown in Fig. 3. The 
effect of kination shows up in the band arHr/2'Kao < f < a^dHkd/'^'^O'a- The spectrum in 
this band is not flat but proportional to yi+(3~2j^)_ Therefore the gravitational waves at high 
frequency are amplifled much more than that in ordinary inflation models. The duration 
of kination depends on the efficiency of reheating. The longer the duration is, the larger 
the amplitude becomes. The gravitational particle production is a least efficient reheating 
mechanism, which may give rise to the amplitude of gravitational waves large enough to 
be detectable. However the energy density of the gravitational waves is constrained by 
the success of the Big Bang nucleosynthesis, fig^ < = 5 x IO'^/iq^ [19i|. Hence the 

efficiency R of the reheating by gravitational particle production has a lower limit. 



R > 0.2{Q^''/5 X 10 
The sensitivity of the advanced LIGO will be Q 



h,Y\H,d/lO'' GeVy 



(42) 



LISA will 



10-'^ at ImHz 



3. 



■gw 



10 
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at 10 ~ lOOHz and that of 



Thus the detectability gravitational waves by these 
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experiments appears unlikely. 

When we calculate the gravitational wave spectrum we assumed that v is constant in the 
inflationary phase. This is because v changes rapidly at the end of inflation in our toy model. 
However, if we consider a model that v changes mildly, we must consider corresponding 
effects on the spectrum. If v is larger than 3/2 in the smaller scales than the COBE bound 
region (fig^ < 7 x 10""(/7o//)^/io ^ IQ-^^Hz < / < IQ-^^Hz H), the slope in the bands of 
the kination becomes looser and the energy density of the gravitational waves may be high 
enough in the detectable bands by LIGO and LISA. 

Instant preheating is so efficient that the duration of kination is much shorter than that 
in gravitational reheating. This results in suppression of the enhancement of the spectrum of 
the gravitational waves at ~ lOOMHz. Hence we find that the spectrum of the gravitational 
waves is quite sensitive to the reheating process after inflation. From Fig. 3 we may obtain 
information regarding the thermal history of the universe after inflation from the spectrum 
of the gravitational waves at ~100MHz. 



V. CONCLUSION 



We investigated the possibility of quintessential inflation by an exponential-type potential. 
In this inflation model, the present accelerated expansion is explained by the same scalar 
field that caused inflation in the past. We considered two extremes of reheating process: the 
gravitational particle production which is least efficient and the instant preheating which 
is most efficient and has not been considered in the context of quintessential inflation. We 
calculated the spectrum of gravitational waves produced during inflation. We found that 
the spectrum of the gravitational waves is quite sensitive to the reheating process. Less 
efficient reheating results in larger amplitude of the spectrum. We emphasize that these 
results are not limited to quintessential inflation and insensitive to the detailed model of 
inflation since the spectrum is essentially determined by the expansion rate during inflation 
and that at the horizon-crossing time. On very large scales corresponding to the Hubble rate 
at an early stage of inflation, it is severely constrained by CMB observations, while on very 
small scales, the dominant constraint comes from BBN which is less stringent. Therefore, 
while it is definitely important to obtain the shape of the spectrum at high frequencies (~ 
100 Hz) or at very high frequencies (~100 MHz), even the detection or nondetection of the 
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gravitational waves would provide us with useful information of the early universe. 

In particular, for detectors with arm length of 1 m, for example, the typical frequency 
range would be 100 MHz. Such tabletop gravitational wave detectors may be constructed 
with much less cost [2^. The implications of such high frequency detectors on the early 
universe would be enormous. 

Acknowledgments 

The authors would like to thank the Yukawa Institute for Theoretical Physics at Ky- 
oto University, where this work was initiated during the workshop on gravitational waves 
(YITP-W-01-16). This work was supported in part by a Grant-in- Aid for Scientific Research 
(No. 15740152) from the Japan Society for the Promotion of Science and by a Grant-in- Aid 
for Scientific Research on Priority Areas (Nos. 14047212 and 14047214) from the Ministry 
of Education, Science, Sports and Culture, Japan. 



[1] P. J. E. Peebles and A. Vilenkin, Phys. Rev. D 59, 063505 (1999). 

[2] M. Peloso and F. Rosati, JHEP. 12, 26 (1999); G. Huey and J. E. Lidsey, Phys. Lett. B 514, 
217 (2001); A. S. Majumdar, Phys. Rev. D 64, 083503 (2001); V. Sahni, M. Sami and T. 
Souradeep, Phys. Rev. D 65, 023518 (2002); K. Dimopoulos and J. W. F. Valle, Astropart. 
Phys. 18 , 287 (2002). 

[3] M. Giovannini, Phys. Rev. D 59, 123511 (1999). 

[4] M. Giovannini, Class. Quantum. Grav. 16, 2905 (1999). 

[5] D. Babusci and M. Giovannini, Phys. Rev. D 60, 083511 (1999). 

[6] A. Riazuelo and J.-P. Uzan, Phys. Rev. D 62, 083506 (2000). 

[7] G. N. Felder, L. Kofman, and A. D. Linde, Phys. Rev. D 59, 123523 (1999); Phys. Rev. D 
60, 103505 (1999). 

[8] L. Kofman, A. D. Linde, and A. A. Strarobinsky, Phys. Rev. D 56, 3258 (1999). 

[9] F. Lucchin and S. Matarrese, Phys. Rev. D 32, 1316 (1985). 
[10] C. L. Bennett et al., lastro-ph/0302207. 
[11] P. Binetruy, Phys. Rev. D 60, 063502 (1999). 



15 



[12] T. Futamase and K. Maeda, Phys. Rev. D 39, 399 (1989). 

[13] T. Chiba and M. Yamaguchi, Phys. Rev. D 61, 027304 (2000). 

[14] T. Chiba, Phys. Rev. D 60, 083508 (1999). 

[15] L. H. Ford, Phys. Rev. D 35, 2955 (1987). 

[16] T. Damour and A. Vilenkin, Phys. Rev. D 53, 2981 (1996). 

[17] B. Allen, Phys. Rev. D 37, 2078 (1988). 

[18] N. D. Birrell and P. C. W. Davies, QUANTUM FIELDS IN CURVED SPACE, Cambridge, 
(1982) 

[19] M. Maggiore, Phys. Rept. 331, 283 (2000). 

[20] S. Kawamura, private communication. 



